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The age of the universe is obtained in a subset of Cardassian models by using 
WMAP data. Cardassian expansion is a modification to the Friedmann equation that 
ahows the universe to be flat, matter dominated, and accelerating, without a vacuum 



, component. Since this model changes the evolution of the universe, we should not a 

\ priori expect the Cardassian age to be the same as the WMAP Friedmann derived 
result of 13.7 ± 0.2 Gyrs. However, in the subset of Cardassian models we consider, 

• we discover that the age of the universe varies from 13.4-13.8 Gyr over the range 

Q . of parameter space we explore, a result close to that of the standard A Cold Dark 

, Matter model. The Hubble constant h, which may also vary in these models, likewise 

\ varies little from the Friedmann result. 
> ■ 

^ : 

0\ . I. INTRODUCTION 

<^ I Over the past five years, Cosmic Microwave Background (CMB) data have shown the 

^ ■ universe to be flat [l|, |1| and supernova data have indicated that it is accelerating 0| . The 
, matter content of the universe falls well short of the necessary energy density to provide a 
'"^1 flat curvature in the standard cosmological model, the Friedmann model. The most popular 
I ' interpretation of this mismatch is that the "missing" density is assumed to be present in 
2 ■ the form of a vacuum (or dark) energy that provides a pressure leading to the acceleration 
to . of the universe. 

\ Alternatively, general relativity may need to^e modified. Modification of the standard 

. ^ ' cosmological model, as in the Cardassian model [5] , allows both the flatness and acceleration 
^ ■ to be accounted for solely due to the matter content of the universe. This modification 
changes the evolution of the universe. Hence the CMB anisotropy spectrum will imply a 
different age of the universe jaL 

The Cardassian model 0, H, 1^ modifies the Friedman equation by adjusting the 

right hand side to be a more general function of the energy density. This function returns 
to the usual Friedmann equation for the early history of the universe, so that ordinary 
nucleosynthesis takes place. However, in the recent past, beginning at a redshift 0(1), this 
function drives the universe to accelerate. Such modifications to the Friedmann equation 
may arise if our universe exists on a brane embedded in a higher dimensional universe jlol |. 
Other proposed modifications to the Friedmann equations include (lll |. 

In section m we will review the Friedmann model and describe how to determine the age 
of the universe. In section UTTl Wilkinson Microwave Anisotropy Probe (WMAP) results are 
reviewed and the choice of appropriate parameters is discussed. In section HVl we will expand 
these techniques to the Cardassian model. We note that, since the customary parameters h, 
Vth, and VLm are not orthogonal and indeed have values that depend upon the cosmological 
model itself, more suitable parameters must be extracted from the data in order to obtain 
the universe age in any model. And finally, results for the Cardassian model are discussed in 
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section In this paper, we obtain the age of the universe only for the Modified Polytropic 
(MP) Cardassian models defined below. 

II. AGE IN THE FRIEDMANN MODEL 

In standard cosmology, the evolution of the universe is governed by the Friedmann equa- 
tion: 

3 ^ ' 

where H is the Hubble parameter and we have dropped the curvature term as experimental 
results are consistent with a flat universe [H . At the curent epoch the critical density is 

= = 1.054 X lO-^/i^ GeV/cm^ (2) 

where subscript refers to the present day and Hq = lOO/i km/s/Mpc. Here, Q = y is 
the ratio of energy density to the critical density, with Qi representing the ratio due to 
the component pt of the density. In the standard picture, an additional component beyond 
matter and radiation is assumed to reach the critical density. This component is taken to be 
a vacuum energy: a cosmological constant A or a time dependent vacuum energy, or scalar 
field known as "quintessence" , that evolves dynamically with time 3 . 



The Hubble parameter can be related to its present day value by: 

H = HoEf{z) (3) 

EpizY = Qoril + Zf + fion^(l + zf + ^]ox(l + zf^^^^^^ (4) 

where f^or, ^om, and VLqx are the current contributions from radiation, matter, and vacuum 
respectively, and wx = Vx/ Px is the equation of state. 

To determine the age of the universe, we integrate Eqn. (jT]). 

, dz 1 

JO [1 + z)Ef[z) 



III. WMAP PARAMETERS 

The WMAP Collaboration has determined values for several parameters by mapping the 
cosmic microwave background relevant ones are listed in Table HI They have determined 
the age of the universe in a A-CDM (cold dark matter) model which assumes that the dark 
energy is a cosmological constant, 

to = 13.7 ± 0.2 Gyr (A-CDM model), (6) 

where the uncertainty is determined by statistical analysis of their Monte Carlo Markov 
Chain results [13]. 

Using the WMAP results for h and f2om in Eqn. ^ yields the correct value, but over- 
estimates the uncertainty. This is due to the fact that the parameters h and fiom are not 
orthogonal; that is, their uncertainties are correlated 0,01- Instead, the set of (nearly) 
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WMAP Parameters 



Hubble constant h = 0.7ll[!;[S;^ 

baryon density VLq^ = 0.044 ± 0.004 

ojb = 0.0224 ± 0.0009 

matter density ^om = 0.27 =b 0.04 

, , _ n 1 Qc.+o-oos 

'^m — L'.±OtJ_QQgg 

acoustic scale Ia = 301 ± 1 

redshift of decoupling Zdec = 1089 ± 1 

TABLE I: Cosmological parameters determined by the WMAP Collaboration from fits to multiple 
experiments Q]. 

orthogonal parameters are uj^ = ^omh"^, = ^obh"^ (where b is for baryon), and the acoustic 
scale Ea- 

Oscillations in the photon-baryon fluid during the early universe led to peaks in the CMB 
power spectrum. The first peak in this spectrum is related to the angle 9a subtended by 
the conformal distance s a sound wave travelled from the big bang until decoupling at the 
surface of last scattering, which has a conformal distance denoted by D. In a flat universe. 



D 



The distances s and D are given by: 



Ho Jo Ef{z 
where Cg is the speed of sound in the photon-baryon fluid: 

1 



1 c 

s = — dz (8) 

-no J^dec ^F{Z) 
1 r^dec 1 

D = irL d^TTT-, (9) 



3(1 + 1^ 



(10) 



An in-depth discussion of a photon-baryon fluid and the CMB may be found in Refs. 0, 



The three parameters oOm, and (.a have nearly orthogonal effects on the CMB power 
spectrum; Ca is related to the position of the first peak, while the other two are related to 
the peak heights. 

The redshift of photon-baryon decoupling z^ec is determined mainly from the overall 
temperature of the CMB, with a WMAP value of 1089 ± 1. The redshift of matter-radiation 
equality z^q is given by: 

5464 ( T \' 

^^^^ 1 -fp,/p^ 0.135 l2.725Ky ^ > 

with a neutrino to photon density ratio of Pi^/p^ = 0.6851 and an overall CMB temperature 
of T = 2.725 K EF 



4 



While several of the above equations may be greatly simplified in the standard Friedmann 
model (with closed forms available for some of the integrals) we will not make these 
simplifications here. Instead, it will be useful to use numerical techniques that can be applied 
below to the Cardassian Model, where the equivalent integrals do not have closed forms. 
From Eqn. we can take h to be an implicit function of Um-, and ^A'■ 

h = h{uJm,uJbjA) (12) 

That is, given observed values for Um, ujb, and £a (from WMAP), we can use numerical 
routines to find the value of h necessary to satisfy = where s and D (dependent upon 
ujrn, (^6, and h only) are given by Eqns. (jH)) & ©, respectively. The value of h thus obtained 
may then be inserted into the time integral, Eqn. (0) (using i^om = ^^m/h'^), to obtain an 
age that is dependent on the three orthogonal parameters. 

This procedure correctly reproduces both the Hubble constant h and age of the universe 
to (including uncertainties) obtained through the WMAP MCMC, listed in Table El Repro- 
duction of the correct uncertainties affirms that the aforementioned parameters are indeed 
nearly orthogonal and will be the appropriate parameters to use in more general models. 



IV. AGE IN CARDASSIAN MODELS 

Cardassian expansion was proposed as a model in which matter alone is sufficient to 
drive acceleration of the universe; in this model, there is no vacuum energy. Instead, the 
Friedmann equation is modified to give a general function of the energy density on the right 
hand side: 

= g{p) (13) 

This function returns to the usual Friedmann equation for the early history of the universe, so 
that ordinary nucleosynthesis and evolution results. However, in the recent past, beginning 
at a redshift 0(1), this function drives the universe to accelerate. 
We see the critical density pc is then defined by the relation: 

9{Pc) = Hi (14) 

where we note the critical density is not, in general, the same as in the original Friedmann 
model (Eqn. (j2))). It will be useful to relate the new critical density and other parameters 
to the old critical density. To prevent confusion, we will use pc to denote the Cardassian 
critical density, while pc will continue to denote the Friedmann critical density. Likewise, 
Vti = y will continue to be defined in terms of the Friedmann critical density. 



In this section, we will look at a subset of the Cardassian models |5l lal known as Modified 
Polytropic Cardassian (MP Cardassian): 



r2 ^T^Gp 



H' = 9{p) 



p \ g(n-l)' 



^ Pcard f 



1/9 



(15) 



where Pcard-, <? > 1, and n <2/'i are free parameters. One of these parameters is fixed by the 
observed value of VLmo-, and a constraint on the remaining two is obtained by matching super- 
nova data (the redshift at which the second term becomes important). The normalization 
is fixed to match the Friedmann model at early times. 
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We will require n < 2/3, so that the second term in the brackets of Eqn. ()15|) increases 
over time, eventually dominates, and then leads to acceleration of the universal expansion. 
In this paper we consider only n > 0; in fact it could be negative as well. The case g = 1 is 
equivalent to the original power law form of the Cardassian model , and the case g = 1 & 
n = is equivalent to A-CDM with the second term representing the vacuum component. 

Taking the geometry to be fiat, the critical density condition Eqn. ()14|1 gives: 



Pc = Pc 



Po 

Pcard 



q{n-l)' 



(16) 



Here the critical density can be smaller than in the original Friedmann case. We will choose 
the parameters so that matter alone is sufficient to provide today's critical density, without 
any vacuum contribution, pom ~ Po = Pc so that 



^Om = 



POm 
pc 



pc 
Pc 



(17) 



The numerical value of Qom will be roughly 1/3 when matching this model to the CMB data 
(see the discussion in the results section below). The condition in Eqn. (fTTj) is satisfied for 
parameters satisfying 

(18) 



pcard — Po 



9(1-") 



or, equivalent ly. 



1 + z, 



card 



^'Om -L 



39(l-n) 



(19) 



After this redshift, the second term on the RHS of Eqn. (|15p begins to dominate and (for 
n < 2/3) the universal expansion will accelerate (even without a vacuum component). 
With the previous constraints, we may rewrite Eqn. (|15|) as: 



H = HoEc{z,q,n) 

Ec{zf = nor{l + zY + ^lom{l + zf + ^loxfx{z) 



where (treating fior as negligible): 



^ox 
fx{z) 



1 ~ ^Om 

{l + zf 
- 1 



l + (Oo^-l)(l + 2;; 



|3cj(n-l) 



1/9 



1 



(20) 
(21) 



(22) 
(23) 



and we now have only two free parameters, q and n. 

Note that at high redshifts z, Ec{z,q,n) Ep{z) for all allowable choices of q and n, 
so that at early times, Ec{z, q, n) is equivalent to the Friedmann case (compare Eqn. (jH)) 
to Eqn. (j^n|) ). Consequently, the evolution of the universe before and around the period of 
decoupling {z > 1000) is the same in the Cardassian model as in the Friedmann model so 
that ordinary nucleosynthesis and oscillations in the photon-baryon fiuid are unaffected. 

The age of the universe may be derived in the same manner as in the standard Friedmann 
model, and we find 

to = Ho-' (24) 
Jo {1 + z) Ec{z,q,n) 
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n 


h 


tn (Gvr) 


1 





0.71 ± 0.04 


13.7 ± 0.2 


1.5 


0.2 


0.72 ± 0.04 


13.6 ±0.2 


2 


0.3 


0.72 ±0.03 


13.6 ±0.2 


10 


0.4 


0.76 ±0.03 


13.4 ±0.2 


100 


0.4 


0.77 ±0.03 


13.4 ±0.2 



TABLE II: Values for the Hubble constant h and age of the universe to for selected SNe la allowed 
parameters in the Cardassian model 0]. 

This integral does not in general have a closed form and must be evaluated numerically. 

Note the CMB peak amplitudes depend upon the parameters Um and Uf,, but are essen- 
tially independent of the late time expansion. Thus, the WMAP fits to these two parameters 
are valid- the Cardassian model does not affect their determination from these peak am- 
plitudes. The values for fiom and Qqi, are not fixed, however, as h is not independently 
determined from the amplitudes. Instead, h may be determined from the peak positions, 
expressed by Ia, which also depends upon the expansion model. In our analysis, we may fix 
LJm, ujb, and Ea to the WMAP observed values and proceed as before. 

Since the sound horizon s in Eqn. (jS)) is evaluated at 2; > 1000, its value is unaffected 
by Cardassian modifications which are important only at late times. Hence s is essentially 
independent of q and n; additionally it has only a minor dependence on h. For a given Ub 
and ujrn, s varies by less than 0.1% over the entire allowed parameter space. However, the 
distance to the surface of last scattering D is evaluated over a period when the Cardassian 
modifications are important so that D has a non-trivial dependence on q and n. 

Since the CMB peaks are generated prior to and during the time of last scattering, 
when the usual Friedmann approximation holds, Ia remains the same, independent of the 
Cardassian parameters q and n. Hence Ia is fixed by the WMAP data, having the same 
value in the Cardassian model as in the standard Friedmann case. Due to Eqn. ((7j), the fact 
that £a and s are unchanged relative to the standard A model constrains the value of D to 
be unchanged as well. However, since D depends non-trivially on q, n, and h, as described 
in the previous paragraph, we must allow all three parameters to vary. In particular, the 
value of h must be allowed to vary in such a way that D remains constant; i.e., 

h = h{um, uJb, ^A, q, n) (25) 

That is, for WMAP observed values for u^, ^^b, and Ea, and for chosen values of the pa- 
rameters q and n, numerical routines can be used to find the value of h necessary to satisfy 
£^ = The chosen parameters q and n and the numerical solution for h may then be 
used to determine the age of the universe using Eqn. (j^ . 

V. RESULTS 

We have obtained the age of the universe appropriate to the MP Cardassian modification 
to the Friedmann equation given in Eq. (fT3j) . Our results are displayed in Figures ^ and |21 
Since the best fit value of the Hubble constant now depends on the new parameters n and q, 
we present a contour plot of h over this parameter space in Figure [T] A contour plot of the 



7 



0.8 



0.6 



0.4 



0.2 








10 

q 



100 



FIG. 1: Hubble constant h contours (in units of 100 km/s/Mpc) for various values of q and 
n in the MP Cardassian model (defined in Eq.(|15|)). based upon central values of WMAP data. 
Uncertainties associated with several cases may be found in Table ^ The_solid lines show the 
supernovae la allowed region of parameter space determined by Wang et al. 



age is shown in Figure |21 The contours were generated using the central values for ojm, ujb, 
and observed by WMAP; uncertainties for several choices of the parameters are shown 
in Table im 

Constraints from SNe la data on the free parameters q and n, which were obtained by 
Wang et al. 0|, are illustrated in both the h and age contour plots. These constraints 
were derived using fixed f2om = 0.3; recall, however, that we have fixed Um in our analysis, 
in which case f2om varies. It would likewise be more appropriate to fix for deriving 
SNe la constraints. The constraints in this case, though, would not significantly differ. For 
n ^ and ~ 1, the MP Cardassian model is nearly equivalent to the A-CDM model and 
^om = 0.3 is valid; hence the SNe la constraints in this region should remain valid. For 
g 3> 1, the luminosity distances to the supernovae are essentially independent of Qom and 
the constraints shown in this region are also valid |3^] . With the fixed f2om constraints valid 
for q ^ 1 and q ^ 1, the intermediate q constraints are unlikely to be significantly different 
from the fixed Um values (the constraints would otherwise have to exhibit fairly pathological 
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FIG. 2: Age of the universe contours (in Gyrs) for various values of q and n in the MP Cardassian 
model (defined in Eq.(|15|)). based upon central values of WMAP data. Uncertainties and SNe la 
allowed region are as described in Figure ^ 



behavior). So while the SNe la constraints shown here are not rigorously correct for our 
analysis, they are fairly close to the "true" constraints and should still be valid. 

Other comparisons of the MP Cardassian model to observational data have also been 
done B H HI H H m, 113. In addition to the SNe la data, Wang tt al. determined 
constraints on the q-n parameter space from the shift in the CMB angular spectrum j^; 
these constraints are weaker than those of the SNe la data, except at g < 2, where the 
upper bound on n is somewhat reduced (see Figure 2 of their paper). Amarzguioui et al. 
have combined CMB and SNe la measurements to determine constraints on the parameter 
space very similar to that of the Wang et al. SNe la constraints, but with a preference for 
lower values of q (see Figure 3 of their paper; the la contour gives q <4). Note, by 
including the CMB spectrum with the SNe la observations, the concerns over the validity of 
the Wang et al. limits discussed in the previous paragraph do not apply to the limits given 
by Amarzguioui et al. 

Values for h and to are shown in Table |n] for several cases of q and n lying within the 
SNe la allowed region. The case g = 1 & n = 0, which is equivalent to the A-CDM model, 
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reproduces the WMAP results with a Hubble parameter h of 0.71 ± 0.04 and an age of 
13.7 ±0.2 Gyr. 

Even over the entire parameter space we considered for the MP Cardassian model (with 
q ranging from 1-100), the age varies very little, ranging from 13.2 to 14.6 Gyr. The SNe la 
allowed band roughly follows the contours and contains values rangin g fr om 13.4 to 13.8 
Gyr. This range is consistent with other constraints on the age |28l. I29L 13(1 . l3l[ l3^ : in 
particular, it is above the minimum of to > 12.6^14 Gyr determined from globular cluster 
ages by Krauss & Chaboyer [3^ as well as a minimum of to > 12.5 ±3 Gyr determined from 
radioisotope studies by Cayrel et al. |3^ . 

While the CMB determination of h depends on the overall model, other techniques mea- 
sure its value more directly. The Hubble Space Telescope (HST) Key Project analyzed 
Cepheids over distances of 60 to 400 Mpc to determine a value of 0.72 ± 0.08 [34]. Mea- 
surements over such short distances are not significantly affected by the variations in the 
cosmological model. While the value of h obtained in Figure 1 does vary significantly over 
the entire allowed (n < 2/3) parameter space, ranging from 0.50 to 0.97, we can see that 
the SNe la allowed band roughly follows the contours and contains values between 0.66 and 
0.78. This allowed band is thus consistent with the HST results. Thus we find the MP 
Cardassian model to be in general agreement with other independent determinations and 
limits for the Hubble constant and age of the universe. 

We now have estimates of the matter density in these models as well. From the value of 
uOm measured by the CMB (given in Table IHl , together with the estimates of the Hubble 
constant given in Figure Q we can extract VL^. For the full n < 2/3 parameter space 
discussed in the last paragraph, we find 0.14 < Vtom < 0.54. The range is far more restricted 
if we consider only those values consistent with the SNe la allowed band: 0.22 < VLom < 0.31; 
this result is in agreement with that of Amarzguioui et al. ji^] iP'Om ~ 0.3 in their analysis). 

In summary, we have shown how to determine the age of the universe in a MP Cardassian 
model from CMB data, which requires the appropriate choice of CMB derived parameters. 
We find that the CMB implied age is entirely consistent with other observational constraints 
for essentially the entire MP Cardassian parameter space examined here (alternatively, ob- 
servational constraints on the age do not constrain the Cardassian parameter space). When 
restricting the parameter space to that consistent with SNe la measurements [8|, the age of 
the universe falls somewhere between 13.4 and 13.8 Gyr, a range that includes the age in a 
A-CDM model. 
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